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O ■ Abstract 
<N 

We propose and investigate a procedure to measure, at least in principle, a positive quantum 

D , 

version of the local kinetic energy density. This procedure is based, under certain idealized limits, 
on the detection rate of photons emitted by moving atoms which are excited by a localized laser 
beam. The same type of experiment, but in different limits, can also provide other non positive- 
0^ , definite versions of the kinetic energy density. A connection with quantum arrival time distributions 
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is discussed. 
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I. INTRODUCTION 



To obtain an expression for the local density of a quantum observable not diagonal in 
coordinate representation, one may look for guidance to the corresponding classical case. For 
a classical dynamical variable, A(q,p), in position- momentum phase space its local density, 
cxa(x), is simply obtained by 

a A (x) = J dpp(x,p)A(x,p) (1) 
= J dqdpp(q,p)5(x — q)A(q,p) 

where p(q,p) is the phase space density. To quantize this expression one can use 

5(x — x) = \x) (x\ (2) 

and consider, for a point x, the operator A(x) = A\x)(x\, or rather one of its many sym- 
metrizations, as a quantum density for the observable A. For a given state \ip), the expec- 
tation value (ip\A(x)\ijj) would then be a candidate for the value of the local density at the 
point x of the observable A. If A is not diagonal in coordinate representation so that it does 
not commute with |x)(a;| , there are infinitely many "combinations" (orderings) to construct 
a quantum density,^ for example, 

I 1/2 |x)(x|l 1/2 (3) 
-(A\x)(x\ + \x)(x\A) (4) 

- (I 1/2 |a;)(x|I 1/2 ) + J (A\x)(x\ + \x)(x\A) (5) 

The noncommutativity of two observables does not mean that there is only one "true" sym- 
metrization of their product. Different symmetrizations may have a perfectly respectful 
status as physically observable and measurable quantities, and different orderings may be 
associated with latent properties that may be realized via different experimental measure- 
ment procedures. They may also be related more indirectly to observables and yet carry 
valuable physical information. An example of this non-uniqueness due to different sym- 
metrizations is the arrival time of a quantum particle at a particular position. Classically 
the distribution of arrival times would be the flux for particles moving in one direction, but 
quantum mechanically different quantizations have been proposed.^ 
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Another important example of this quantum non-uniqueness for a single classical quantity 
is the kinetic energy density^^ in which case A becomes the kinetic energy operator, 
A = p 2 /2m. There is no unique definition of a quantum "kinetic energy density", in spite 
of the relevance of the concept in several fields. The Thomas-Fermi theory provides an 
early example of a possible realization and application. In density functional theory, it 
enters as one of the terms of the energy functional to determine the electronic structure of 
atcas, .decu.es, 8oli ds 0r fennionic gases, see e,. Refs. ffl In thia context it haa been 
used in particular to define a local temperature and identify molecular sites most reactive 
to electrophilic reagents.- The kinetic energy density also plays a key role in partitioning 
molecular systems into fragments with well defined energies and virial relations ,^ 10 ' 11 or to 
define "intrinsic shapes" of reactants, transition states and products along the course of a 
chemical reaction.— It is moreover a basic quantity in quantum hydrodynamic equations.^*^ 

In all these applications, different quantum versions of the kinetic energy density have 
been used. The most commonly found cases are the three quantizations considered above, 
or suitable generalizations thereof. They satisfy different properties and there have been 
intensive discussions which one is best, but clearly they all are useful in different ways 
and for different purposes. However, not so much emphasis has been placed on possible 
procedures of how to measure them. 

Many of these arguments and controversies already can be seen in the simple case of the 
kinetic energy density of a free particle in one dimension. Real, Muga and BrouarcU^ studied 
and compared three versions of an operator, t(x), for the kinetic energy density at a point 
x, associated with different quantizations, namely 

\ p5(x-x)p 

T {X)= 2m ' (6) 



p 2 p 2 

— S(x — x) + 5{x — x) 

2m 2m 



(7) 

Gill 

?W(x)+T®{xj\ . (8) 



The second operator follows from the quantization rule of Rivier.^ The corresponding den- 
sity {f^ 2 ')t is given by its, generally time dependent, expectation value and may in principle 
be obtained operationally by a weak measurement of the kinetic energy post-selected at 
position x . 16 i 17 i 18 The third one, which is the average of (f^')t and (f^ 2 ^, corresponds to 
Weyl's quantization rule. An indirect way to measure (f^)t for free motion was described 
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by Johansen,>l£ who noticed that the second time derivative of the expectation value of \x— x\ 
is proportional to (t^)<. 

In this paper we will provide an operational interpretation of the first expression which, 
incidentally, is the only positive one among the three, and in fact among a much broader 
family of quantizations.— We will use for this end a simple model, originally devised to study 
time of arrival measurements.— The basic physical idea is to send atoms in their ground state 
through a region illuminated by a perpendicular laser beam and to measure the resulting 
fluorescence rate of photons. 



II. KINETIC ENERGY DENSITY, FLUORESCENCE AND ATOMIC ABSORP- 
TION RATE IN AN IMAGINARY POTENTIAL BARRIER 

The description of photon fluorescence from moving atoms, which are excited by a local- 
ized laser beam, is based on the quantum jump approac h 21 ' 22 ' 23 and has been discussed in 
detail elsewhere . 20 ' 24 ' 25 We will only summarize the results which are relevant for the present 
investigation and assume a "lambda" configuration of three atomic levels in which the laser 
couples levels 1 and 2 with Rabi frequency Q, whereas level 2 decays with inverse life time 7 
predominantly and irreversibly to a ground (sink) state 3j2£ For a laser on resonance with the 
atomic transition, atoms for which no photon is detected ( "undetected atoms" ) are governed 
by the following effective Hamiltonian 

H = f/2m + hl{x) {|2)(1| + |1)(2|} - ^ 7 |2)(2| , (9) 

where Q(x) is the position dependent Rabi frequency. The evolution of the wave function for 
undetected atoms simplifies to a one channel Schrodinger equation if h'j is large compared 
to the kinetic energy, with a complex imaginary potential 27 

u(x) = -ih^. (io) 

In this "low saturation regime" the undetected atoms, whose number is proportional to 
the norm-squared of the wave-function, are in the ground state most of the time, and 
the temporal probability density for a photon detection ("detection rate") is given by the 
decrease of the number of undetected atoms, 

n(i) = -d(m\m)/dt = -j[w*)MW(*)> ■ (11) 
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This is the basic operational quantity obtained in the modeled experiment. We will relate 
it, or its normalized version 

n ~ (t) = nm (12) 

to ideal quantities for the freely moving atom, unperturbed by the laser, by taking certain 
limits. 

First, we consider a square laser beam profile. Then the imaginary potential becomes a 
barrier of height 

V = htt 2 /2^ (13) 
located between x = and x = L, and the effective Hamiltonian for undetected atoms reads 

where X\o,l] is 1 inside the laser illuminated region and zero outside. The absorption, or 
detection, rate is found from Eq. and is given by 

n(t) = — y o dx\^{x,t)\ 2 - (is) 

To obtain the time development of a wave packet coming in from the left, we solve first the 
stationary equation 

H<P k = E k <j) k . (16) 

Using standard matching conditions, the energy eigenfunctions <pk in the barrier region 
< x < L for a plane wave coming in from the left with momentum hk are given by 

(A + {k)e iqx + A_{k)e~ iqx ) (17) 

a / /TV \ ' 
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with q 2 = k + 2imV/h and 

A + (k) = Kq±k)e^ 
±y ' 2kqcos{qL)-t{k 2 + q 2 )sm{qL)' 1 ' 

Writing the initial state as a superposition of eigenfunctions with positive momenta, we 

obtain 

r°° ~ 9 
ij)(x,t)= I dki){k)e- %hkt/2m (t) k {x). (19) 
Jo 

Inserting this into Eq. ()15|) yields the absorption rate. It is also useful to define the auxiliary 
freely moving wave packet 

^ / (ac,t) = -=/ dkip{k)e- ihkt/2m e %kx . (20) 



2tt 



We will now relate the absorption rate to an ideal distribution by going to the limit V — > oo. 
When V is increased, more and more atoms are reflected without being detected, but nor- 
malizing the result a finite distribution is obtained in the limit, even though the absorption 
probability vanishes eventually due to total reflection. For large V, V ^> h 2 k 2 /2m, one has 
in leading order 



q ~ y2imV/h 2 , 
2k 



Q 

2k 
Q 



,(,., ~ -^=— (e igx + e iq{2L - x) ) , < x < L. (21) 
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Integrating over x and neglecting the terms which vanish exponentially, the absorption rate 
becomes 

n(t)~ —= I dk dk'^*{k)ij{k')e ih{k2 - k ' 2)t/2m kk'. (22) 

nmyrnV Jo Jo 

This expression is independent of the barrier length L as a result of the large V limit, so the 
same result is obtained with an imaginary step potential —iVQ(x) or with a very narrow 
barrier. 

The normalization constant is given by JU(t)dt ~ 2hk (mV)- 1/2 , where k = 
J \ip(k)\ 2 k dk, and the normalized absorption rate is 

fi roo poo _ „ „ 

n N (t)~—L L / dk dk'^*(kWk')e in{k k )t/2m kk' . (23) 
2nmk Jo Jo 

With the freely moving wave packet ipf{x ) t) this can finally be rewritten in the form 

n "(*) - ^(Mt)\kS(m^f(t)) • (24) 
vnKo 

Now, the right hand side is just the expectation value at time t of the kinetic energy density 
evaluated at the origin! Thus, with pq = ko% the initial average momentum we have 
obtained, in the limit V — ► oo, 

lim n N (t) = -(¥»(x = 0)) t . (25) 

V^oo p 

Note that the averages are computed with the freely moving wave function and that one 
obtains the kinetic energy density at an arbitrary point a by shifting the laser region, i.e. 
by replacing [0, L] by [a,L + a] in Eq. (fTf]) . 



Remark: Instead of normalizing the absorption rate of Eq. (|15|1 by dividing by a constant 
one can normalize it on the level of operators,— which preserves it as a bilinear form. 
However, in this case the result depends on the constant to which one normalizes. It has 



been shown in Ref. 



25l that if the constant is chosen as 1 then operator normalization of Eq. 



(|15[) leads for V — > oo to the arrival time distribution of Kijowski.— Now, since the time 
integral of (f^ (x = Q)) t equals po/2 = hko/2, it is suggestive to choose po/2 as normalization 
constant. Following the approach of Ref. |30 operator normalization of n(t) then leads to 
(7^(1 = 0))t in the limit V — > oo. 



III. KINETIC ENERGY DENSITY FROM FIRST PHOTON MEASUREMENT 
AND DECONVOLUTION 

In the preceding section we had V = Q 2 /2 , ~y, with much larger than the kinetic energy. 
Therefore, the limit V — > oo implies a simultaneous change of Q and 7. Experimentally, the 
Rabi frequency Q is easy to adjust, but not the decay rate 7. To overcome this problem we 
therefore describe a procedure that allows to keep the value of 7 fixed. 

We again consider the moving-atom laser model but now for the limit Q — > 00 and 
7 = const. In that case, the simplified description of the evolution of the wave function by 
means of the imaginary potential U(x) is not feasible, and one has to solve the full two- 



channel problem 
done in Refs. 2G 



or the three- level atom with the Hamiltonian in Eq. ©. This has been 



25l . and normalizing with a constant the resulting photon detection rate 



Ilw(t) becomes 



U N (t) ~ — — / dkdk'r(k)ip(k') e ^ k2 - k ' 2 W 2m ^ , ow • (*3) 

w 27rmk Jo V 1 y ' 7 + ih(k 2 - k' 2 ) / 'm v ; 



For h'j large compared to the kinetic energy of the incident atom, Eq. (J23|) is recovered, 
but for finite 7 there is a delay in the detection rate. This can be eliminated by means of a 
deconvolution with the first-photon distribution W(t) for an atom at rest . 2 ^ 2 ^ The ansatz 

n(t) = U id (t) * W{t) (27) 

yields in terms of Fourier transforms 

^ = Wi\ ^ 
W(u) 
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with 20 i 31 



W{v) \ fi 7/ fi 7 fi2 

~ 1 + , ft -> oo. 29 

7 



Inserting this and the Fourier transform of Eq. (Jzbj) into Eq. (|28|) . the resulting ideal 
distribution, after performing the inverse Fourier transform, reads 

7i OO POO 

U id {t)~——- dk dkY{k)^{k')e ih{k2 ~ k ' 2)t/2m kk', (30) 



27imk Jo 

which is the same expression as the absorption rate of Eq. ([23)1 . obtained here operationally 
for fixed 7. Naturally, 



n^)~-(f«(x = o)) t (31) 



holds as before. 



IV. CONNECTION WITH QUANTUM ARRIVAL TIME DISTRIBUTIONS 

Here we briefly discuss a formal connection between the kinetic energy density r^(x) 
given in Eq. (|24|) and Eq. (|25|) and the arrival-time distribution of Kijowski 29 , 

U K (t) = ^ f (t)\k^5(x)k^\^ , (32) 

at x = 0. For wave packets peaked around some ko in momentum space, the operator k 1 ! 2 
acting on tpf in Eq. (J3*2|) can be expanded in terms of (k — ko), 

k 1/2 = k l Q /2 + h 1/2 (k - k ) - \k^'\k - k ) 2 + 0((k- k ) 3 ) . (33) 

In the following we take k to be the first moment of the momentum distribution, k = 
J \ip(k)\ 2 kdk. Inserting the expansion in Eq. ()33|) into Eq. (|32|) yields in zeroth order a very 
simple result, 

U K (t)=v \M^t)\ 2 + O(k-k ) , (34) 

i.e. the particle density times the average velocity vq = koh/m. To first order in (k — ko) 
one obtains the flux at x = 0, 

U K (t) = J(0,t) + O((k-k ) 2 ) , (35) 
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where 

J(P,t) = 7^(Mt)\$m + 6(x)k)\1> f (t)) , (36) 
and to second order the expression 

U K {t) = J(0, t) + J-A(0, t) + 0{(% - k f) (37) 

where 

A(0,f) = (fW(x = 0)) t - (r«(x = 0)) t . (38) 

For states with positive momentum, which we are considering here, the first order, namely 
the flux, is correctly normalized to one, and so is the second order since the time integral 
over A is easily shown to vanish. This difference only provides a local-in-time correction to 
J that averages out globally. Its quantum nature can be further appreciated by the more 
explicit expression 

i = * n« (39) 

2p Hmp ox 

which shows the inverse dependence on mass and momentum and the explicit quadratic 
dependence on %. 



V. DISCUSSION 



In Fig. 1, operational and ideal kinetic energy densities are compared for a coherent su- 
perposition of two Gaussian wave packets with different mean momenta. They are prepared 
in such a way that their centers of mass arrive simultaneously at the origin. This enhances 
the interference among different momentum components and the differences between the 
distributions. As seen in the figure, the differences between various versions of the quantum 
kinetic energy density may be quite significant. While (r^{x))t is always positive, (r^ 2 '(x))t 
can become negative in classically forbidden regions for stationary eigenstates of the Hamil- 
tonian, a fact that has been used by Tachibana 12 to define molecular and reaction shapes. 
It is perhaps less obvious that this quantity can also be negative as a result of free motion 
dynamics, as seen in the figure. 

The main contribution of this paper has been to point out that, under certain idealiza- 
tions and limiting conditions, fluorescence experiments can lead to an operational approach 
to kinetic energy densities. A positive kinetic energy density can be obtained from the "mea- 
sured" signal in a strong laser limit. With A the difference of the positive density 
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and the density (t 2 ) 4 of Rivier, an interesting relation was found for A with the ideal time- 
of-arrival distribution of Kijowski and with the flux. The latter two can also be obtained 
operationally under appropriate limits. In a recent review by Ayers, Parr and Nagy 8 - on 
the kinetic energy density, one of the suggestions for future research was the need to study 
and understand this quantity A better. In a completely different context from the present 
work, and suitably generalized to three dimensions, A plays a major role in Bader's theory 
to separate a molecular system into meaning ful fragments.-^ii Note that, if A = 0, (r^) t 
becomes also equal to the other two densities considered. Therefore this condition implies a 
certain "classicality" or coalescence of the multiple quantum possibilities. If A or its integral 
over some volume are zero, a fragment can be defined with a well defined kinetic energy and 
virial relations. It is quite striking that, in the second order expansion of the arrival time 
distribution of Kijowski, A = implies that LT^ becomes the flux, which is, as we know, the 
quantity that plays the role of a arrival time distribution in classical mechanics. 

In summary, we think that these results clarify the status, as physically meaningful 
physical quantities, of several versions of the local kinetic energy densities, and may stimulate 
experimental research on quantum kinetic energy densities and as well as on arrival times. 
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FIGURE CAPTION (Fig. 1) 

Comparison of kinetic energy densities at x = 0: (f w ) (solid), (f( 2 )) (dashed), (f( 3) ) 
(dotted) and the operational quantity p IlN(t)/2, see Eqs. (fT2|) and for V = 1.9 hfJST 1 , 
L = 0.21 fim (triangles) and V = 950.0 h /is" 1 , L = 0.42 fim (circles). The initial wave packet 
is a coherent combination ip = 2~ 1 ^ 2 (ipi + ip 2 ) of two Gaussian states for the center-of-mass 
motion of a single caesium atom that become separately minimal uncertainty packets (with 
Ax i = Ax2 = 0.031/im, and average velocities (v)i = 18.96 cm/s, (v) 2 = 5.34 cm/s at 
x = and t = 2 /is). The mass is 2.2 x 10 25 kg and p = 2.67 x 10 _26 kg m/s. 
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Muga, Seidel and Hegerfeldt, Figure 1 
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